Scaling laws of turbulence and heating of fast solar wind: the role of density 

fluctuations 
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Incompressible and isotropic magnetohydrodynamic turbulence in plasmas can be described by an 
exact relation for the energy flux through the scales. This Yaglom-like scaling law has been recently 
observed in the solar wind above the solar poles observed by the Ulysses spacecraft, where the 
turbulence is in an Alfvenic state. An analogous phenomenological scaling law, suitably modified to 
take into account compressible fluctuations, is observed more frequently in the same dataset. Large 
scale density fluctuations, despite their low amplitude, play thus a crucial role in the basic scaling 
properties of turbulence. The turbulent cascade rate in the compressive case can moreover supply 
the energy dissipation needed to account for the local heating of the non-adiabatic solar wind. 
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The interplanetary space is permeated by the solar 
wind pj], a magnetized, supersonic flow of charged parti- 
cles originating in the high solar atmosphere and blowing 
away from the sun. Low frequency fluctuations of solar 
wind variables are often described in the framework of 
fully developed hydromagnetic (MHD) turbulence 0,3- 
The large range of scales involved, spanning from 1 AU 
(~ 1.5 X 10® km) down to a few kilometers, make the so- 
lar wind the largest "laboratory" where MHD turbulence 
can be investigated using measurements collected in situ 
by instruments onboard spacecraft Q. MHD turbulence 
is often investigated through the Elsasser variables = 
V ± (47rp)~^/^ b, computed from the local plasma veloc- 
ity V and magnetic field b, p being the plasma mass den- 
sity. In terms of such variables, MHD equations can be 
rewritten as dtz^ -\- • \7z^ — —VP/p + diss, where 
P is the total hydromagnetic pressure, and diss indi- 
cates dissipative terms involving the viscosity and the 
magnetic difFusivity. As in the Navier- Stokes equations 
for neutral fluids, the nonlinear terms z^ • Vz* cause 
the turbulent energy transfer between different scales, at 
high Reynolds numbers where dissipative terms can be 
neglected. However, in the MHD case, they couple the 
two Elsasser variables, so that the Alfvenic MHD fluctu- 
ations z*, propagating along the background magnetic 
field, are advected by fiuctuations z^ propagating in the 
opposite direction. The presence of strong correlations 
(or anti-correlations) between velocity and magnetic fiuc- 
tuations, along with a nearly constant magnetic inten- 
sity and low amplitude density fiuctuations, is usually 
referred to as Alfvenic state of turbulence, and implies 
that one of the two modes z* should be negligible, mak- 
ing the nonlinear term of MHD equations vanish for pure 
Alfvenic fiuctuations. In that case, the turbulent energy 
transfer should also disappear ,4]. Alfvenic turbulence 
is observed almost ubiquitously in fast wind. This holds 
both in the ecliptic fast streams, and in the high latitude 
wind blowing directly from the sun coronal holes (1, @ . 



As pointed out in [1], the observation of Alfvenic state 
turbulence in the solar wind represents therefore a para- 
doxical "contradiction in terms" . 

MHD turbulence however satisfies an important ana- 
lytical relation, which is the equivalent for magnetized 
fiuids of the Kolmogorov or the Yaglom relations. Un- 
der suitable hypotheses, it has been shown that the 
pseudo-energy fiuxes {£) through the scale i of the in- 
crements of the Elsasser fields Az^{£) = z^{x-\-£)—z^{x) 
follow a linear scaling relation 



(1) 



Here, represents the component of the incre- 

ment Az^ along the direction £, and are the dissipa- 
tion rates per unit mass of the pseudo-energies {\z^\'^)/2 
((•) indicates space averages). The scaling law ^ has 
been recently observed experimentally in polar wind Q 
and in the ecliptic plane [3, [13] • The confirmation of the 
scaling law ([T]) is an important step towards the solution 
of the apparent paradox of the Alfvenic turbulent state, 
because it unambiguously shows that an MHD cascade 
is present, maybe with a weak transfer rate, despite the 
strong velocity-magnetic fields correlations. 

Relation ([TJ , which is of general validity within MHD 
turbulence ^s not always realized in the solar wind ob- 
servations Q. Indeed, it is possible that local character- 
istics of the solar wind plasma do not always satisfy the 
assumptions required for ^ to be valid, namely large- 
scale homogeneity, isotropy, and incompressibility. The 
role of anisotropy of solar wind turbulence, which is ex- 
pected to be important, is not considered in this letter. 
Density fiuctuations in solar wind have a low amplitude, 
so that nearly incompressible MHD framework is usu- 
ally considered However, compressible fluctua- 
tions are observed, typically convected structures charac- 
terized by anti-correlation between kinetic pressure and 
magnetic pressure [l^ . Properties and interaction of the 
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basic MHD modes in the compressive case have also been 
considered in the past EH- In the present paper, 
we show that density fluctuations, despite their very low 
amplitude, play a central role in the turbulent energy 
transfer, and that a phenomcnological scaling law ob- 
tained by taking into account density fluctuations is ob- 
served in a much larger proportion of fast solar wind. 
We also show that the turbulent dissipation can account 
for a large fraction of the local heating causing a slower 
than expected decrease of temperature with distance. A 
flrst attempt to include density fluctuations in the frame- 
work of fluid turbulence was due to Lighthill [l^. He 
pointed out that in a compressible energy cascade, the 
mean energy transfer rate per unit volume ev pv^ /£ 
should be constant in a statistical sense (v being the 
characteristic velocity fluctuations at the scale £), ob- 
taining V ~ {1/ pY^'^ . Fluctuations of a density- weighted 
velocity field u = p^/'^w should thus follow the usual 
Kolmogorov scaling i. The same phenomcnolog- 

ical conjecture can be introduced in MHD turbulence 
by considering the pseudo-energy dissipation rates per 
unit volume Cy = pe^, and introducing density- weighted 





lo' 


m 


10 






[km 


10 




10^ 








10^^ 






% 


10^ 




10> 



1996, day 23 




10" 



10° lO' 
T [hours] 



10^ 



FIG. 1: One example of the mixed third order compressible 
pseudo-energy flux (r) as computed from the Ulysses data 
during days 23 to 32 of 1996. The incompressible flux Y^{t) 
in the same time window and a hnear fit are also indicated. In 
this case, both compressible and incompressible fluxes obey a 
Yaglom-like law. 



Elsasser fields, defined as 
of the Yaglom-type relation 



W^{1) = (|Aw±|"Au;^) {p) 



should then hold for the density-weighted incre- 
ments Aw*(^). Note that we have defined the 
flux W^{t) so that it reduces to Y^{t) in the case of 
constant density, allowing for comparisons between the 
compressible scaling ^ and the purely incompressible 
one ([T]). Despite its simple phenomcnological deriva- 
tion, the introduction of the density fluctuations in 
the Yaglom-type scaling ([2]) seems to describe correctly 
the turbulent cascade for compressible fluid (or mag- 
netofluid) turbulence. The law for the velocity field has 
been observed in recent numerical simulations fl^, [l^ . 

We will now study the cascade properties of compres- 
sive MHD turbulence from solar wind data collected by 
spacecraft Ulysses. In order to avoid as far as possi- 
ble variations due to solar activity, or other ecliptic dis- 
turbances such as slow wind sources, coronal mass ejec- 
tion, current sheets, we concentrate our analysis on pure 
Alfvenic state turbulence observed in high speed polar 
wind. We use here measurements from Ulysses space- 
craft in the first six months of 1996. Such period was 
characterized by low solar activity, so that solar origin 
disturbances werebv almost absent. Moreover, the space- 
craft orbit was at high and slowly decreasing heliolati- 
tude, from about 55° to 30°, and presented small vari- 
ations of the heliocentric distance r, from 3 AU to 4 AU. 
Since the mean wind speed (v) in the spacecraft frame is 
much larger than the typical velocity fluctuations, and is 
nearly aligned with the radial direction R, space scales I 
can be viewed as time scales r, related through the Tay- 
lor hypothesis by ^ = — {vji) t. We then used 8 minutes 
averaged time series of both Elsasser variables {t) and 



The equivalent density p{t) = rip + inuc (obtained as the sum of proton 
density and 4 times He density), to compute the density- 
^ weighted time series w^{t). From this time series we cal- 

f ^ £ (2) culate the increments Aw^ (r) = {t + t) — {t) for 

3 different time lags r, and the third-order mixed structure 

functions W^It) = (|Aw='=(r)pAw^(T))t by time aver- 
aging {■)t over windows of fixed duration t. The same pro- 
cedure has been also been used to calculate the quantities 
Y^{t) using the time series of the Elsasser fields z^{t). 
In order to eliminate instationarities, heliolatitude and 
heliocentric distance changes, and to explore the wind 
properties locally, averages are computed over a moving 
window of about 11 days, consisting of 2048 data points. 
Accuracy of third order moments estimate 19] was tested 
with such sample size [2d| . We found that the third-order 
structure functions W^{t) computed from the Ulysses 
data show a linear scaling 



(3) 



during a considerable fraction of the period under study. 
In particular, we observed linear scaling of W~^{t) in 
about half of the signal, while W^{t) displays scaling 
on about a quarter of the sample. As comparison, the 
corresponding incompressive scaling law for Y^{t) was 
only observed in a third of the whole period, considerably 
smaller than the compressible case Q. The portions of 
wind where the scaling is present are distributed in the 
whole period, and their extensions span from 6 hours up 
to 10 days. The linear scaling law generally extends on 
about 2 decades, from a few minutes to one day or more. 
For the compressible scaling, the two fluxes W^{t) co- 
exist in a large number of cases. This does not hold for 
the incompressive scaling, where in general the scaling 
periods for the two fluxes Y^{t) are disjoint. 

Figure [1] shows one example of both mixed third-order 
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FIG. 2; Top panel: an example of the third order compress- 
ible pseudo-energy flux W^{t) during days 1 to 10 of 1996. 
Bottom panel: W~(t) for days 66 to 75 of the same year. In 
both panels, the corresponding incompressible fluxes F*(r) 
(no scaling present) and a linear fit are displayed. 



structure functions W^(t) and Y^(t) computed in the 
same 11 days windows where the scaling was observed. 
Figure [5] shows two more examples of scaling, observed 
both for W^(t) and W^~(r), in two different time win- 
dows. The W^{t) scaling extends over 2 decades, while 
W~(t) behaves linearly on the whole range of scales con- 
sidered here (3 decades). In the last example, the scal- 
ing is not present for the incompressible fluxes V^ij). 
This example shows that the inclusion of compressible ef- 
fect through the density-weigthed fluctuations improves 
the scaling ([2]) and modify the energy cascade. The 
scaling relation © also allows a direct estimate of the 
pseudo-energy transfer rates in the compressible case. 
A fit of the linear law ([3]) provides the local values of 
the amount of pseudo-energy transferred from large to 
small scales by the turbulent MHD cascade. This was al- 
ready measured in the incompressive case [2l| , so that 
it is possible to compare the transfer rates in the two 
cascades. The mean values, computed over the 46 ob- 
served scaling cases at different radial distances from the 
sun, (± their dispersion, in [J kg^^ sec^^]) for the com- 
pressible cascade are e+ — 3668 ± 1900 (29 cases) and 
e~ = 3536 ± 2500 (17 cases). Both values are consider- 
ably larger than the corresponding values for the incom- 
pressive case (e^ = 182 ±73, 24 cases, and ej = 156 ±50, 
11 cases [2l[). This result shows again that the cascade 



in the solar wind is strongly enhanced by density fluctu- 
ations, despite their small amplitude. Note that the new 
variables are built by coupling the Elsasser fields with 
the density, before computing the scale dependent incre- 
ments. Moreover, the third order moments are very sen- 
sitive to intense field fluctuations (intermittency), that 
could arise when density fluctuations are correlated with 
velocity and magnetic field. Similar results, but with 
considerably smaller effect, were found in numerical sim- 
ulations of compressive MHD |22|] . We should point out 
that experimental values of energy transfer rate in the in- 
compressive case had been also estimated with different 
tecniques from different datasets [1, [l3| . Those values are 
not in agreement with the present (incompressive case) 
results. However, the different nature of wind (ecliptic vs 
polar, fast vs slow, at different radial distances from the 
sun) makes such comparison only indicative. 

An interesting open question is the problem of the so- 
lar wind heating. The first models of solar wind assumed 
an adiabatic cooling due to spherical expansion of plasma 
blowing out of the sun. This would result in a radial de- 
crease of the proton temperature T(r) ^ with ^ = 
4/3. On the contrary, spacecraft measurements [13] have 
shown that the temperature decay is slower than the adi- 
abatic prescription, with ^ £ [0.7,1]. This implies that 
some local heating mechanism is present. One stand- 
ing hypothesis is that the heating could be provided by 
energy dissipation occurring at the small scales of a tur- 
bulent cascade [1,11^113]. Using equation (P), the rates 
at which the incompressible turbulent pseudo-energy is 
transported down the scales, and eventually dissipated at 
small scale, can be measured directly from data. This has 
recently be used to investigate whether or not a turbulent 
cascade can heat the solar wind. Results were however 
not conclusive. In fact, the measured incompressible dis- 
sipation rate of pseudo-energies can only account for up 
to 50 % of the solar wind heating [l^, Hlj . 

Figure [3] shows the radial profiles of the pseudo-energy 
transfer rates for both the compressive and incompres- 
sive cascades. In the same figure, we show the profiles 
of the heating rates needed to obtain the observed tem- 
peratures, as estimated from heating models ^ |2l|, 
and from the measured temperatures (the two different 
values refer to the different estimates of the tempera- 
ture obtained from Ulysses instruments). It is evident 
that, while the incompressive cascade cannot provide all 
the energy needed to heat the wind, the density fluc- 
tuations coupled with magnetohydrodynamic turbulence 
can supply the amount of energy required. This evidence 
shows the importance of the density fluctuations in po- 
lar, fast solar wind turbulence, confirming that it should 
be considered as an example of compressive fully devel- 
oped MHD turbulence. Note that, since in a few sam- 
ples we measured both e+ and (T in the same period, 
the values of the energy e = (e+ + e~)/2 and cross- 
helicity tn = (e^ — e )/2 transfer rates can be disen- 
tangled. From the values obtained, it is clear that the 
cross-helicity contribution, indicating the importance of 
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FIG. 3: Radial profile of the pseudo-energy transfer rates 
obtained from the turbulent cascade rate through the Yaglom 
relation, for both the compressive and incompressive case. 
The solid lines represent the radial profiles of the heating rate 
required to obtain the observed temperature profile. 

the Alfvenic state of turbulence, can vary from a neg- 
ligible fraction (less than 1 %) to a considerable 25 % of 
the energy contribution. Since its amplitude does not ap- 
pear to be correlated with the observation of the cascade, 
Alfvenicity seems not to play a crucial role in the cascade 
at the observed scales. This would be in agreement with 
previous analysis of solar wind turbulence anisotropy, 
where the Alfvenic contribution to the field fluctuations 



is small [25|.[26|. 

To summarize, we used the density- weighted Elsasser 
fields to show for the first time that a phenomenologi- 
cal compressive Yaglom-like relation is verified to a large 
extent within the solar wind turbulence. This implies 
that low amplitude density fluctuations play a crucial role 
for scaling laws of solar wind turbulence [22| . This obser- 
vation also confirm the recent results for the Kolmogorov 
4/5-law from numerical simulations of compressible tur- 
bulence while no experimental evidences from real 
fluids had been found so far. This could be attributed 
to the incompressible nature of flows in ordinary fluids 
accessible to laboratory experiments. Here in fact, we 
present the first experimental observation of relation ([2]) 
in real systems. Using solar wind data, we have access 
to a sample of weakly compressible MHD turbulence in 
nature. Scaling law is found to be quite common and ex- 
tends on a large range of scales, indicating not only that 
a nonlinear MHD cascade for pseudo-energies is active 
in the solar wind turbulence, but also that compressible 
effects are an important ingredient of the cascade. We 
point out that the observed departures from the scal- 
ing law could be due to presence of inhomogeneity and 
anisotropy in the solar wind [28| . The compressive cor- 
rections to the cascade also cause the transfer of a consid- 
erably larger amount of energy toward the small scales, 
where it can be dissipated to heat the plasma locally. The 
role of anisotropy in the solar wind turbulent cascade still 
remains an open question. 
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